A SYSTEM EXHIBITING TOROIDAL ORDER 



A. B. Harris 

Department of Physics and Astronomy, University of Pennsylvania, Philadelphia PA 19104 

(Dated: June 29, 2010) 

A two dimensional system of discs upon which a triangle of spins are mounted is shown to undergo 
a sequence of interesting phase transitions as the temperature is lowered. We are mainly concerned 
with the 'solid' phase in which bond orientational order but not positional order is long ranged. As 
the temperature is lowered in the 'solid' phase, the first phase transition involving the orientation or 
toroidal charge of the discs is into a 'gauge toroid' phase in which the product of a magnetic toroidal 
parameter and an orientation variable (for the discs) orders but due to a local gauge symmetry these 
variables themselves do not individually order. Finally, in the lowest temperature phase the gauge 
symmetry is broken and toroidal order and orientational order both develop. In the 'gauge toroidal' 
phase time reversal invariance is broken and in the lowest temperature phase inversion symmetry 
is also broken. In none of these phases is there long range order in any Fourier component of the 
average spin. A definition of the toroidal magnetic moment T; of the ith plaquette is proposed such 
that the magnetostatic interaction between plaquettes i and j is proportional to TiTj. Symmetry 
considerations are used to construct the magnetoelectric free energy and thereby to deduce which 
coefficients of the linear magnetoelectric tensor are allowed to be nonzero. In none of the phases 
does symmetry permit a spontaneous polarization. 

PACS numbers: 75.25.+z,75.10.-b,75.50.Ee,77.80.-e 



I. INTRODUCTION 

The theoretical analysis of toroidal ordering in 
electric [ll-Q and magnetic [l||3Il| systems has recently 
been investigated. Examples of such states in which 
plaquettes of spins assume a chiral configuration have 
been known for some time.Q However, it is widely be- 
lieved that toroidal magnetic order should always be sub- 
servient to the primary order parameter, a Fourier com- 
ponent of the average spin. Here we address the possibil- 
ity of defining a toroidal order parameter for a system in 
which it is the primary magnetic order parameter. The 
major problem is to identify a situation in which there 
is ferrotoroidicity but there is no nonzero Fourier com- 
ponent of average spin. In this paper we consider a two 
dimensional system of toroidal spin plaquettes, such as 
those shown in Fig. [1] which exhibits the desired behav- 
ior. 
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that it is aligned either parallel or antiparallel to its local 
axis fixed in the plane of the disc, as shown in Fig. [5J The 
intraplaquette dipolar interactions are strong enough so 
that at temperatures of interest the spins in each plaque- 
tte come to thermal equilibrium in one of the two degen- 
erate ground states as shown in Fig. [1] The magnetic 
dipole moment of each spin is mimicked by a pair of oppo- 
site charges (iQ), as shown in Fig. [2] If r is the distance 
of the charges from the center of the plaquette, then the 
magnitude of the dipole moment is p = 2Qrsinx. The 
orientation of the disc is defined by the angle (f) between 
the x-direction and the dashed line fixed on the plaque- 
tte. 
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FIG. 1: Toroidal plaquettes with 2, 3, or 4 moments. 
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II. THE MODEL 



FIG. 2: A plaquette with three spins in one of their dipolar 
ground states. Each magnetic dipole is represented by a pair 
of charges ±Q whose positions are fixed by the angle x- 



We consider a system of microscopic circular discs 
(confined to lie in the x-y plane) which contain three 
spins in a triangular configuration as in the center panel 
of Fig. [T] Each spin has a large single ion anisotropy so 



We now calculate the interaction energy Vab of pla- 
quette A whose center is at the origin and plaquette B 
whose center is at R ee {X,Y) = (i?cos^, i?sin\E'). To 
calculate Vab we invoke the multipole expansion of the 
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electrostatic energy (which is proportional to the desired 
magnetic dipole-dipole energy) between the charges on 
the two plaquettes. This energy is a sum over terms of 
the form (l/i?)F„((/)A, 0b, 5')z", where z = r/R. The 
positions are functions of rexp(±ix) and the symmetry 
of the plaquette indicates that the result can only involve 
factors of r'^ exp(±3ix). Also this term must be odd in 
both XA and xb- Accordingly, the lowest order nonzero 
term is of order / R. So at this order 

Vab = A((/.A,0s,*)QAsin(3xA)Qi3sin(3xs)'^V^^ 

(1) 

where Qa [Qb) is the amplitude of the charges on pla- 
quette A (B). We only need this for small x, so 



Vab = 9A{(l)A,(l)B,'f)QAXAQBXBryR\ 



mAmBr'^A{(f)A,(l)B,^) 



(2) 



where rrii — 2Qirxi is the ith magnetic moment. For 
the present situation it seems reasonable to define the 
toroidicity of the ith plaquette such that 



so that 



Vab = A{cbA,cbB,'^)TATB/R'' , (3) 



= (3/2)TO,r2 = 3Q^x^r^ ■ (4) 



Note that this definition, unlike the standard one,[7| has 
the property that the magnetostatic interaction energy 
between two toroids is proportional to product of their 
toroidal 'strengths' T^. To get A(0a, '/>b, ^) we may con- 
sider the lowest order terms in the expansion of Vab hi 
powers of XA and xb- Thus 



A{(t>A,(t>B,'^) 



R' 



9QaQb dxAdxB (6!)dr6 
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The phase space of this model is specified as follows. 
Each plaquette is characterized by its center of mass po- 
sition R inside a two dimensional box in the x-y plane 
with respect to which the plaquette has mirror symmetry. 
The ith plaquette can assume an orientation specified by 

and has a toroidal strength Ti which is an Ising-like 
variable: Ti = ±(3/2)|mi|r'^, according to which of the 
two ground states the spins occupy. We assume an ori- 
entationally independent interaction between plaquettes 
which for concreteness we take to be the Lennard- Jones 
potential, 
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where a and e are constants. Q Because the system is two 
dimensional, the plaquettes can not develop long range 
order characteristic of a three dimensional solid, pj In- 
stead the system can develop various behaviors interme- 
diate between a conventional solid and an isotropic liquid. 



As shown in Ref. 1^ the system exhibits a 'solid' phase 
in which bond-orientational order is truly long range but 
position correlations exhibit power-law decay. The 'solid' 
melts either directly or indirectly (via a hexatic phase) 
into an isotropic liquid phase. Here we focus on the tran- 
sition as the temperature is lowered through the value T/ 
at which the 'solid' phase appears and assume that kTj is 
much larger than the magnetic interactions between pla- 
quettes. Even in the 'solid' phase, there is no true long 
range positional order. Then the spin correlation func- 
tion will involve an average over correlations between a 
spin at the origin and spins in a distant plaquette. Since 
the average spin in a plaquette is zero, the spin correla- 
tion function, like the positional correlation function, can 
not display long range order. In contrast, because there 
is long-range bond orientational order, lowering the tem- 
perature can lead to phase transitions due to the interpla- 
quette dipolar interactions. It is the purpose of this paper 
to analyze the symmetry of the resulting ordered phases. 
In a likely scenario we find that the 'solid' undergoes two 
such phase transitions. At the first (higher temperature) 
transition we find that time reversal symmetry is bro- 
ken and at the second spatial inversion symmetry is also 
broken. Since the underlying system does not have long- 
range positional ordering, the spin correlation function 
itself is never long ranged. To substantiate this picture 
it is necessary to analyze the interplaquette interactions 
and thereby verify that they lead to Ising-like transitions. 



A. Interplaquette Interaction 

To calculate the interaction between two plaquettes, 
each of which is confined to the x-y plane, we assume (see 
Fig. [2) that plaquette A has charges <jaQa, where a a = 
±1, at (positions relative to the center of the plaquette) 



r{nA,crA) = [a:;(nA, cr^), j/(nA, cr^)] 



(7) 



where 



x{nA,cFA) = rcos(cr^XA + 2nA7r/3 + '/'a) 
= [r cos(xa) cos(2nA7r/3 + 0a) 

-CTA?'sin(xA)sin(2nyi7r/3-|-0A), (8) 

y{nA, <7a) = r sin{aAXA + 2nA7r/3 + (Pa)] 
= aAr sin(xA) cos(2nA7r/3 + 0a) 

+r cos(xa) sin(2nA7r/3 -1- 0a)] (9) 

for HA = 1,2,3 and similarly for plaquette B. 

Then the interaction energy between the two plaque- 
ttes is 
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= 1 (TAO"B=±l 



R' 



-2i?^ [x{nA,<yA) - x{nB,crB)] 
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+2Rr] [y{nA, a a) - y{nB,(7B)] 
-2r(nA, (Ta) ■ r{nB,crB) + 2/ 



-I -1/2 



, (10) 



where ^ = cos(*) = X/R, and 77 = sin* = Y/R. We 
drop the last term in since it can not contribute to 
a term involving sin(3xA) sin(3xB)- We want to expand 
this in powers of z and now we work to first order in XA 
and xb by setting sin(xx) = Xx and cos(xx) = 1, where 
X is either A ot B. Thus 



Vab 
where 



QaQi 
R 



^ ^ a^c7B7^-l/^ (11) 



TZ = I + 2^z cos(2n^7r/3 + 4>a) — cos(2nB7r/3 + 4>b) 



OAXA sin(2nA7r/3 + (?!)yi) + (JbXb sin(2n_B7r/3 + ^s) 

2r]z i^AXA cos(2nA7r/3 + 0a) 

obXb cos(2nB7r/3 + (/^b) 

+ sin(2nA7r/3 + 0a) — sin(2nB7r/3 + 0b) 



2z^ cos(2nA7r/3 + 0a) — ctaXA sin(2nA7r/3 + 0a) 



X cos(2nB7r/3 + 0s) — gbXb sin(2nB7r/3 + 0_b) 



2z^ fTAXA cos(2nA7r/3 + 0a) + sin(2nAvr/3 + 0a) 



X \^bXb cos(2nB7r/3 + 0b) + sin(2nB7r/3 + 0b)J (12) 

We now analyze how Vab depends on i. e. how 
it depends on ^ and t]. Note that ^ and 77 each carry a 
factor of z. However, factors that do not depend on 77 
and f carry factors of 1 or z^, so that can be no terms 
at order which are cubic in 77 or ^. Therefore Vab can 
only depend on through the argument G^*. In addition, 
for Vab to be invariant under a global rotation, it can 
only depend on differences in angle. Also the result can 
only be a function of 30a and 30b because it must be 
invariant under 0^ 0^ + 27r/3. Thus A must be of the 
form 

A(0A,0s,*) =a + &cos(6^'-30A-30B) 
+ccos(30A — 30b) + dsin(6* — 30a — 30b) 
+esin(30A - 30b) ■ (13) 

One can view the system looking either along the positive 
z axis or along the negative z axis. Comparing these two 



views, one sees that a positive rotation is equivalent to 
changing the sign of the charges followed by a negative 
rotation. But the interaction energy is invariant under 
charge conjugation. So 



A(0A,0B,4') 



a + 6 cos(6\I' — 30A — 30b) 
+ccos(30A — 30b) • 



(14) 



Furthermore, consider what happens if we average the 
interaction over the orientation 0^ of one of the plaque- 
ttes. This superposition of charges leads to a uniformly 
charge-neutral plaquette. This argument tells us that a 
in Eq. (fT4| must be zero. To identify the coefficients h 
and c it suffices to evaluate 



A(0,0, *) = 6cos(6*) + c . 



(15) 



The explicit evaluation of A(0,0, 4') is carried out in 
the Appendix. Then use of Eq. implies the result 



A = i^^cos(6vE'-30A-30B) 
--^cos(30A -30b) . 



(16) 



It is important to note a local symmetry. With the in- 
teractions so far postulated, the Hamiltonian is invariant 
under the local transformation 



(17) 



This is a nontrivial symmetry that indicates that the two 
configurations shown in Fig. [3] have the same energy 
at leading order in the multipole expansion. Note that 
changing the sign of Qi is equivalent to changing the sign 
of the toroidal moment T;. As a result of this local gauge 
symmetry it follows from Elitzur's theorem^llj that even 
though there is long range order in the variable Q sin(30), 
there is no long range order in either Q or sin(30). 



III. PHASE TRANSITIONS OF THIS MODEL 

Here we give a more detailed analysis of the phase tran- 
sitions within this model. We assume that the magnetic 
anisotropy energy that aligns the spins along a fixed di- 
rection in each plaquette and the isotropic interactions 
of Eq. ^ between plaquettes are dominant. These as- 
sumed interactions do not depend on the orientation of 
either of the interacting plaquettes. (With this assump- 
tion Elitzur's theorem applies.) Accordingly, as the tem- 
perature is lowered, this two dimensional system will un- 
dergo a phase transition at a temperature T/ into a 'solid' 
phase with long-range bond-orientational order, but no 
long range positional order. JT^] It is obvious that in this 
phase both spatial inversion and time reversal symme- 
tries are maintained. 

As one further reduces the temperature, the interpla- 
quette dipolar interactions come into play and can cause 
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of lattice constant a) assumes the form 
q V 



-[cA:r + i.(q)]y(q)y((q)' 



(20) 






FIG. 3: Top: two interacting plaquettes (both with — ir/2 
and Q = +1). Bottom: same as the top configuration ex- 
cept that the right-hand plaquette is rotated through an an- 
gle A(j} — TT and sign of the charge is now Q = —1. These two 
configurations have the same interaction energy truncated at 
order . This is the gauge invariance of Eq. (|17|) . 



order to develop consistent with the local gauge symme- 
try. To see what sort of order develops we introduce the 
appropriate gauge invariant variables 



X, = I^,cos(30,), = T,sin(30,) 



(18) 



Because bond orientational order is maintained, we can 
treat each molecule as being surrounded by a hexagon 
of neighboring plaquettes and the orientation of this 
hexagon of neighbors is maintained over the entire sys- 
tem. Accordingly, we can define 4* as being measured 
relative to the direction between the central plaquette 
and one of its neighbors. So we may take 6^'/(27r) to 
be an integer for all nearest neighbor interactions. We 
do not consider further neighbor interactions in view of 
how rapidly the intcrplaquette interaction falls off with 
separation. For simplicity we work as if we have a two di- 
mensional triangular lattice. Thus we analyze the model 
with orientationally-dependent interactions 



Vab = 



10, hotaTb 



cos{3(f)A) cos(30b) 



io,62or^ra • /o , ^ • /o^ . 

sm(30A) sm(3(/>i3; 



(19) 



So we have a two dimensional anisotropic rotor model 
which is in the same universality class as the two di- 
mensional Ising model. Such models have been widely 
studied, [l^fisj To analyze the phase transitions within 
this model we invoke mean-field theory, within which the 
Landau free energy, in terms of the Fourier transforms 
of Xi and Yi (temporarily assuming a triangular lattice 



at quadratic order, where c is a constant of order unity, 
and /i(q) and i^(q) are the Fourier transforms of the po- 
tential: 

/i(q) = 2A[cos(aqj;) -I- 2 cos(a(7j;/2) cos(\/3a(7j,/2) ] 
iy{q) — 2B[co8{aqx) + 2cos{aqx/2) cos{\/3aqy/2) [2,1) 



where 



A = 



10, 170T2 
32i?7 



B = -- 



10,620T2 
32i?7 



(22) 



where T — (3/2)mr^. Thus as the temperature is lowered 
the system will develop long range "ferro" order (order at 
zero wave vector) in the variable Y = Tsin(3(/>), but not, 
as was said, in either T or sin(3(/)) separately. This type 
of order is illustrated in Fig. U) This transition occurs at 
a temperature of order kTjj = 6\B\ w 2000rV-R'^- 






FIG. 4: Long range order with unbroken gauge symmetry. 
Note that the top-most plaquette is in the gauge transformed 
state according to Eq. (|17|) . Here we illustrate the case when 
the order parameter Q sin(3(?!)) is negative. 



We now discuss whether time reversal (T) symmetry or 
spatial inversion (P) symmetry is broken in this phase. 
As a preliminary, note that the single particle density 
matrix assigns the probabilities p/2, p/2, and 1 — p to 
the states S, S', and S", respectively, where S is the state 
of the top-most plaquette in Fig. 21 S' is the state of 
a plaquette in the bottom row of Fig. SI and S" is the 
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completely disordered state. The interpretation of this 
density matrix is shown in Fig. [5l where we see that 
time reversal symmetry is broken but inversion symmetry 
is maintained. (If the dipole were electric dipoles, then 
inversion symmetry would be broken.) 



1 




1 



state we have antiferrotoroidicity. Here we are mainly 
interested in displaying a ferro state, so we take v to be 
positive. This final ordering transition which breaks lo- 
cal gauge symmetry will occur at a temperature of order 
kTiii — 6w, which we assume to be much smaller than 
kTij. This transition is also in the same universality class 
as the two dimensional Ising model. At this transition 
spatial inversion symmetry is broken. The symmetry of 
the various phases is summarized in Fig. |6l 

new long- < T > ™, • ^ bond 

J < T sin 6> 

range order <- cin rfi ^ ■ ^ j.- 

& < sin 9 > orientation 

symmetry p p j P T 

^ . hexatic or 

, lerro- gauge nonmagnetic . , 

toroid ^ forofd ^ .J,, ^ j^^^ 

Till Ti 

FIG. 6: Symmetry of the various phases. Here the symme- 
tries are T — time reversal symmetry, P — spatial inversion 
symmetry, and 6 = continuous rotational symmetry. 



FIG. 5: The ordered component of the density matrix incor- 
porating states S' and S". This state is not time reversal in- 
variant but is inversion symmetric because spins being pseudo 
vectors do not change their orientations under inversion. 

In this phase (which we call the 'gauge toroid' phase), 
because of the gauge symmetry, there is no long range or- 
der in the orientation of the plaquettes. This is a result 
of our assumption that the interaction between plaque- 
ttes is independent of their orientations. Of course, this 
assumption is not consistent with the three-fold symme- 
try implied by the existence of the spins. Accordingly, 
we now take proper account of this three-fold symme- 
try by introducing small bulges in the plaquettes at the 
locations of the three spins. This will give rise to an 
interplaquette interaction which breaks the local gauge 
symmetry. For simplicity we take this gauge symmetry 
breaking interaction to be of the form 

I^' = ucos(6^'-3(/)a-3<?!)b) . (23) 

Now we again set 6^" / {2tt) to be an integer and we only 
need to consider interactions involving sin(30), so effec- 
tively 

V' = -usin(3</>A)sin(30B) . (24) 

Now what happens depends on the sign of v. If w is pos- 
itive, then we have a ferro arrangement of plaquettes, so 
that all plaquettes are in the same state (either as those 
in the bottom row of Fig. |3]or as that in the top row of 
Fig. S]). Because toroidicity and orientation are strongly 
coupled, this state is ferrotoroidal. If v is negative, then 
we have an antiferro arrangement of plaquettes into the 
so-called 'root-3' structure discussed recently in connec- 
tion with charge ordering in lutetium ferrite.Ti?] In this 



In Ref. @ toroidicity has been discussed in connec- 
tion with the magnetoelectric effect. However, instead 
of using the symmetry of the crystal (see Ref. llQl . they 
used the symmetry of free space to obtain a simplified 
relation between the toroidicity and the linear magneto- 
electric tensor. Here we invoke the symmetry of the two 
dimensional system to write the magnetoelectric free en- 
ergy as a function of the electric field E and the magnetic 
field H as 

•^ME = ^A[Hx:Ey — HyEx] + BHzEz 

+C[H,E, + HyEy]^ , (25) 

where is the z-component of the toroidicity, T. We 
may check that this form is consistent with the symmetry 
of the system, keeping in mind that spin and the mag- 
netic field are both pseudo vectors but T is a real vector. 
Accordingly, both and the factor in the large brackets 
are odd under the mirror z — )• — z. Also all the terms are 
invariant under a rotation about the z axis. Thus 
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= -C . (26) 
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These elements of the hnear magnetoelectric tensor are 
only nonzero in the ferrotoroidal phase. The combina- 
tion of the three-fold axis and the x-y reflection plane 
guarantee that the spontaneous polarization is zero in all 
these phases. (If the discs were asymmetric with respect 
to this mirror, then a spontaneous polarization along z 
would be allowed in the ferrotoroidal phasej^^) 



IV. DISCUSSION AND CONCLUSION 

The nonmagnetic 'solid' to gauge toroid transition at 
Til was analyzed assuming no gauge breaking interac- 
tions. We argue that the introduction of the small gauge 
breaking interaction will not qualitatively modify the 
phase diagram of Fig. [6l because it will take a finite inter- 
action to orientationally order the plaquettes. Similarly, 
including higher order terms in the multipole expansion 
will not alter our conclusions as long asz = r/i?^l. 
One possible difficulty in constructing the system of discs 
analyzed in this paper is that it may be difficult to achieve 
equilibrium statistics within the manifold of the dipolar 
spin ground states. But perhaps it is not crucial that the 
gap between the two spin ground states and the qualita- 
tively different excited spin states is very large. It would, 
of course, be extremely interesting to observe the linear 
magnetoelectric effect in a system such as this. 
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Appendix A: Evaluation of A((;/>a, 0s, 5") 



Write 



Vab 



QaQi 
R 



nA,nB = l CAO-B 
2 



1 + Az + BizaAXA 



+B2Z(JbXb + Ciz aAXA + C2Z (JbXb 
-1-1/2 

+ Dz'^(7A<yBXAXB + Ez^ 



(Al) 



Here 



A = 2^[cos{2nATT/3 + (j)A) - cos{2nBTT/3 + (j)B)] 

+2r] [sin(2n^7r/3 + 4>a) - sin(2nB7r/3 + 0_b)] 
Bi = — 2(_sm{2nA'T/3 + (j)A) + 27] cos{2nATr/3 + (j)A) 
B2 = 2^sin(2nB7r/3 + 0b) — 277C0s(2nB7r/3 -f 0b) 
Ci = 2 sin(2n^7r/3 0a) cos(2nB7r/3 + 0_b) 
— 2cos(2nA7r/3 + 0a) sin(2nB7r/3 + 0b) 



C2 = 2cos(2nA7r/3 -I- 0a) sin(2nB7r/3 -f 0b) 

—2 sin(2nA7r/3 -|- 0a) cos(2nB7r/3 -I- 0b) = — Ci 
D = — 2 sin(2nA7r/3 0a) sin(2nB7''/3 + 0b) 
—2 cos(2nA7r/3 + 0a) cos(2nB7r/3 + 0b) 
E = — 2 cos(2nAvr/3 + 0a) cos(2nB7r/3 + 0b) 

—2 sin(2nA7r/3 + 0a) sin(2nB7r/3 + 4>b) = D . 

(A2) 

Thus, when the sums over cta and ctb are performed, we 
have 



Vai 



^QaQi 
R 
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■ E 



l + Az + Dz' 



"1 "^/^ r 

Dz^XAXB 
2{1 + Az + Dz^) 

3[BiB2z2 + (^^2 - Bi)Ciz^ - Clz^]xAXB 



A{l + Az + Dz^Y 



(A3) 



We have to evaluate this at order z^ . According to Eq. 
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A(0a,0b,^) = ^^ 



z=0 , 



E 



Dz' 



2(1 + Az + i:iz2)3/2 

3[BiB2^' + {.B2 - Bi)Ciz^ - Clz^] 
4(l + ylz + Z3z2)5/2 



(A4) 



Thus 




(-3/2!)i^3 



(-3/2!)L>2A2 



2(-7/2!)(2!) 2(-9/2!)(l!)(2!) 
3(-5/2!)BiB2i:'2 



2(-ll/2!)(0!)(4!) 4(-9/2!)(2!) 
3{-^/2\)BiB2DA^ 3(--5/2!)BiB2^ 
4(-ll/2!)(l!)(2!) 4(-13/2!)(0!)(4!) 
3[-b/2\){B2~Bi)CiDA 

4(-9/2!)(l!)(l!) 
3(-5/2!)(i?2-i3i)CiA3 

4(-ll/2!)(0!)(3!) 
3(-5/2!)Ci2D 3(-5/2!)C2A2 1 



4(-7/2!)(l!) 4(-9/2!)(0!)(2!) 
12 

ri=l 



(A5) 

Bi) separately. 
A(0,0,*) = (4/9)[A6(0,0,*) + AA(0,0,^')] ,(A6) 
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where we number the terms with {B2 
To implement Eq. (1151) we write 



where Ag is the term in A proportional to BiB2A^ and 
AA contains the remaining terms written in Eq. (jASp . 
Note that Ag is the only term which is sixth order in ^ 
and 77 and is therefore the only term which can contribute 
to the term in Eq. (|14p involving G^*. Thus we write 



Ag = a + /3cos(6* - - 3(/)b) 

+cg cos(30A - 3(/)b) (A7) 

and AA is a constant which we can evalute by setting 
= 0. 



AA = a' + c' cos(30A - 30b) 



(A8) 



We first evaluate Ag for (f)A — (t>B — ^ with ^I* arbitrary. 
It is convenient to introduce the notation {x, y), where a: 
and y assume the values 1, c, or s, for unity, the cosine 
function and the sine function. The first argument is 
that for plaquette A and the second is that for plaquette 
B. Thus (c, cs) = cos(2nA7r/3) cos(2nB7r/3) sin(2ns7r/3) 
and (1, c^) = cos^(2nB7r/3). We have the sums over n: 

^(cos(2n7r/3) = ^ sin(2n7r/3) cos'=(2n7r/3) = 

n n 

sin2(2n7r/3) = ^ cos2(2n7r/3) = 3/2 

n 71 

cos3(2n7r/3) = - ^ sin^(2n7r/3) cos(2n7r/3) ==3/4 

n 71 

^cos2(2n7r/3)sin2(2n7r/3) = 3/8 

n 

Y cos'*(2n7r/3) = ^ sin'*(2n7r/3) = 9/8 . (A9) 

n n 

With these understandings (and with the sums over 
UA and ns implied) we write 



(AlO) 



SO that 

56 = -64K(s,l)-77(c,l)][^(l,s)-77(l,c)] 
x(^e[(c,l)-(l,c)]+77[(s,l)-(l,s)]^ 
= -64[e'(s, s) + 7,2(c, c) - c) - ^77(C, s)] 
x(^e[(c,l)-(l,c)]+r;[(s,l)-(l,s)]j 

= -64^2(5, s)]|^,74[_4(^3^^)_4(^^^3)j 
-12r;2^2(^^^)[(^^^)_(^^^)j2j 

-64[,72(c,c)]k^[(c,l)-(l,c)]4 



+6eV[(c,l)-(l,c)]^[(sM) + (l,.2)] 
+771(5, l)-(l,s)]^j 

+64^77(5, c) ^4f77[(c,l)-l,c)]3(s,l) 
+4773^[(c,l)-(l,c)][(s3,l) + 3(s,s2)]j 

+64e77(c, s) ^-4e^/[(c, 1) - (1, c)]3(l, s) 
-4773C[(c,l)-(l,c)][(l,s3) + 3(s^s)]^ 

= 77'(^-64(c,c)[(s,l)-(l,5)f^ 

+77V (^256[(As') + (s',s4)] 
-384(c,c)[(c,l)-(l,c)]2[(s2,l) + (l,s2)] 
+256[(s, c) + (cs)][(c, 1) - (l,c)][(s3, 1) + 3(s- 

+77Y(^768(s^s2)[(c,l)-(l,c)]2 

-64(c,c)[(c,l)~(l,c)r 

+256[(s2,c) + (c,s2)][(c,l)-(l,c)]3j . 

When simplified this is 

56 = r)'^[-64i6){cs^,cs^)]+5Ur]^^\s\s^) 
-76877V[-2(c's2,c2) + (cs2,c3)] 
+512?7V[-(s1 c") + 3(5^0, s^c) - 3(s2, s^^^) 
+76877'C'[(s'c2, s^) - 2{s^c, s^c) + (s^, s^c^)] 
-64?72^1-4(c4, c") + 6(c3, c3) - 4(0^, c^] 
+512772e'[-3(s'c2, c2) + 3(s2c, c3) - d 

= 77^-384] 
-76877^^^^ 



51277*^2 j 



76877^^4 I 



6477^^^ -4-2 
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- 512774^2 
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+ 51277^^"^ -3 



21677*' + rr[1296] - 77^r[1944] , 



(A12) 



This is 



^6 = 1^-108(^2 + 77^)3 

-108(r?^ - 1577*^2 ^ 15^2^4 _ ^e-,^ 

= 108[cos(6*) - 1] . (A13) 
This leads to the result 

93,555 



- COS(30A - 3(t>B) 



128 

+ cos(6* - 3(i)A - 3(/>s) 



.(A14) 



We now evaluate the other sums. But since A A does 
not depend on "Jf, we set ^ = 1 and 77 = 0. Then we have 
(again the sums over ha and ub are left implicit) 

^2)3 = 5i = -8[(s,s) + (c,c)p 
= -8[3(s^c,s^c) + (c^c^)] 



= -18 . (A15) 



If <S2 = D'^A'^, then, with the sums implicit, we have 

^2 = 16[(s^s2) + 2(sc,sc) + (c^c2)] 

x[(c,i)-(i,cr 
+(s^sV) + (c^c4)] 





3' 


2 


■3" 




_~4_ 


+ 


4_ 



32 1^ 
36 . 



(A16) 



If iSs = ^ DA^, then, with the sums implicit, we have 
S3 



32h(..)-(cc)][(c,l)-(l,c)r 
64[2(c^c2)-3(c^c3)+2(c^c4) 

3 



= 64^4 
= 324 . 



(A17) 



If ^4 = J2BiB2D^, then, with the sums implicit, we 
have 



^4 = -16(s,s)[(s,s) + (c, c)] 



-32 



-18 . (A18) 



If 55 = '£BiB2DA'^, then, with the sums implicit, we 

have 

<S5 = 32(s,s)[(s,s) + (c,c)][(c,l)-(l,c)]2 
= 32(s^s2)[(c^l)-2(c,c) + (l,c2)] 
= 32[(s2c2, s^) - 2(cs^ cs^) + s^c^)] 



= 64 



(A19) 



If .Sy = '^B2CiDA, then, with the sums implicit, we 
have 

Sr = -16(l,s)[(s,c)-(c,s)][(c,c) + (s,s)][(c,l)-(l,c)] 
= 16[-(s, sc) + (c, [(c, c) + (s, s)] [(c, 1) - (1, c)] 
= 16[-{s\s^c) + {c^cs^Mcl) - {l,c)] 



(- 


00 


2 


■3' 




3" 




4_ 
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_~4_ 



= -18 . 



(A20) 



J^BiCiDA = Ss = -16{s,l)[{s,c)-{c,s)] 
x[(s,s) + (c,c)][(c,l)-(l,c)] 
= -J2B2C1DA (A21) 

If ^9 = E^2Ci^^, then, with the sums implicit, we have 

59 = 32(l,s)[(s,c)-(c,s)][(c,l)-(l,c)]3 
= -32(c,s2)[(c,l)-(l,c)]3 
= -32[(c^ s^) - 3(c^ cs^) + 3(c^ c^s^)] 



-32 
-162 . 
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74 
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+ 3 


2 
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(A22) 



^ B^CrA' = Sw = mis, 1) - 77(c, 1)] [(s, c) - (c, s)] 
X |^C[(c,l) + (l,c)]+77[(s,l) + (l,s)]j 

= -^S2CiA3 (A23) 

If »Sii = "^CiD, then, with the sums implicit, we have 

5n = -8[(s,c)-(c,s)]2[(c,c) + (,s,,s)] 

= -8[(s^ c^) + (c^ s^)](c, c) -f 16(cs, cs)(s, s) 



= -8 2 



-2 



= 18 . (A24) 



If iSi2 = CiA^, then, with the sums implicit, we have 

5i2 = 16[(s,c)-(c,s)]2[(c,l)-(l,c)]2 

= 16[(s2, + (c2, s2)][(c^ 1) - 2(c, c) + (1, c')] 



9 



32[(s2c2, + {c\ c^) - (s^c, c3) - {c\s^c)] 



32^ 

108 



- 2 



(A25) 



Thus 



V- A n T^ 135 945 25,515 



8 8 
945 



64 



945 



945 25,515 25,515 



8 32 
135 2835 

93,555 



128 



cos(6^') - 



32 



2025 
128 



(A26) 
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